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ABSTRACT 

We present the Combined Array for Research in Millimeter Astronomy (CARMA) ATLAS 3D 
molecular gas imaging survey, a systematic study of the distribution and kinematics of molec- 
ular gas in CO-rich early-type galaxies. Our full sample of 40 galaxies (30 newly mapped 
and 10 taken from the literature) is complete to a 12 CO(1-0) integrated flux of 18.5 Jy km 
s -1 Q, and it represents the largest, best-studied sample of its type to date. A comparison of 
the CO distribution of each galaxy to the g — r color image (representing dust) shows that 
the molecular gas and dust distributions are in good agreement and trace the same underly- 
ing interstellar medium. The galaxies exhibit a variety of CO morphologies, including discs 
(50%), rings (15%), bars+rings (10%), spiral arms (5%), and mildly (12.5%) and strongly 
(7.5%) disrupted morphologies. There appear to be weak trends between galaxy mass and CO 
morphology, whereby the most massive galaxies in the sample tend to have molecular gas in 
a disc morphology. We derive a lower limit to the total accreted molecular gas mass across 
the sample of 2.48 x 10 10 M , or approximately 8.3 x 10 8 M per minor merger within the 
sample, consistent with minor merger stellar mass ratios. 

Key words: galaxies : elliptical and lenticular, cD - galaxies : kinematics and dynamics - 
galaxies : evolution - galaxies : ISM - radio lines : galaxies - galaxies : star formation - ISM 
: molecules - surveys 



1 INTRODUCTION 

The bimodality in the morpholog ical classification of galaxies has 
long been known jHubbld[l926h . Late-type galaxies (LTGs) typ- 
ically have exponential discs, spiral structure and blue colours. 
Early-type galaxies (ETGs) are more ellipsoidal, with smoother 
isophotes and red colours. One scenario for producing this galaxy 
bimodality is via mergers. Simulations have shown that a merger 
between two LTGs o ften creates an ETG (iToomre & Toomrell 19721 ; 
ISnringel et all l2005). ETGs are also much more likely to be found 
on the "red sequence" portion of the color -magnitude diagram 
(CMP: iBaldrv et alJl2004l : iFaber et al.ll2007h . presumably due to 
the loss of most of their star-forming material. 

ETGs have been shown to be deficient in star formation 
relative to LTGs at a given mass dVisvanathan & SandageNl977l ; 
lBoweretal.ll 19921). and therefore should be molecular and atomic 
gas-poor (iLees et alil99lh . Recently, it has however become clear 
that they are not devoid of co l d gas, con t aining reservoirs o f 
dust (e.g. lHawarden et all 1 198 ll: llural 1 19861; Ixnapp et all [l989h . 
neutral atomic gas (e.g. Knappetal. 1985; Sage & Welch 2006; 



Oosterloo et al. boich an d molecular gas (e .g. ISage & Wrobell 
19891 ; IWelch & Sagdl2003l ; ICombes et all 2~~ 



EQQ2). 



The ATLAS 3D sample is a complete, volume-limited sample 
of 260 local ETGs brighter than Mr — —21.5, covering distances 
to 42 Mpc with some restric tions on declination and Galactic lati- 
tude dCappellari et alj2011al hereafter Paper I). The survey was de- 
signed to help us understand ETG formation and evolution. Optical 
integral-field spectroscopy with SAURON on the William Herschel 
Telescope (WHT) has been obtained over the central 41" x 33" re- 
gion of all ATLAS 3D galaxies, revealing their internal stellar kine- 
matics, stellar population properties and ionized gas distributions 
and kinematics. The ATLAS 3D sample has also been completely 
observed in CO J=l-0 and J=2— 1 with the I nstitut de Radioas- 
tronomie Millimetrique (IRAM) 30m telescope (I Young et al.l201lL 



* email: kalatalo@ipac.caltech.edu 
t Dunlap Fellow 

* 18.5 Jy km s _1 is drawn from the IRAM 30m survey, and as the inter- 
ferometric data indicate that the 30m molecular gas masses are generally 
underestimates, the true limit is likely slightly higher. 



hereafter Paper IV), and 65% of the sample has been observed 
in H I with the We s terbork Synthesis Radio Telescope (W SRT; 
iMorganti et al.ll2006l : IOosterloo et al]l2010l : ISerra et al.ll20ll here- 
after Paper XIII). For studies of the cold and warm interstellar 
medium (ISM), the ATLAS 3D sample is thus one of the largest and 
best-observed samples of local ETGs available. 

Of the 260 ETGs in the ATLAS 3D sample, 56 were detected 
in CO (22% detection rate; Paper IV). However, despite emerging 
evidence that a non-negligible fraction of ETGs contain molecular 
gas, little is known about its origin an d evolution, althoug h various 
scenarios have been put forth (see e.g. lDavis et al.ll20lTb[ hereafter 
Paper The three most prominent scenarios are outlined here. 
First, it is possible that the gas survived the galaxy's transforma- 
tion into an ETG and was not completely consumed subsequently 
through secular evolution. Second, the molecul ar gas may have ac- 
cumu lated internally through stellar mass loss jFaber & Gallagher] 
Il976h . In the ATLAS 3D sample, it is expected that each galaxy pro- 
duces on average 0.1 M© yr " 1 of ISM through mass loss from the 
old stellar population alone Third, it is possi- 

ble that the molecular gas was accreted from an external source, 
through tidal stripping events, cold accretion and/or minor mergers 
with gas-rich compani ons. The m i nor merger ( 1:4 - 1:10 mass ra- 
tios) rate calculated bv lLotz et all J2008L [2oTTh suggests that there 
are approximately 0.12 minor mergers per galaxy per Gyr, so the 
ETGs of the ATLAS 3D sample as a whole should have undergone 
a total of « 30 minor mergers in the last Gyr. It is also possible 
that this molecular gas is leftover from major mergers that have 
taken place, but minor mergers dominate in number. Determining 
the origin of the molecular gas therefore requires a detailed anal- 
ysis of its distribution and kinematics (e.g. angular momentum), 
only available from interferometric maps, and a comparison to the 
stellar kinematics (as done in Paper X). 

How the properties of the molecular gas depend on the en- 
vironment of the host galaxy is an open question. For instance, 
galaxies in clusters may well be unable to acquire cold external gas, 
specifically atomic hydrogen (H I), due to the presence of the hot 



1 The current paper presents in detail the observations used in Paper X to 
discuss the origin of the CO. 
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intracluster medium (ICM) and associated ram pressure stripping 
dCavatte etal.l lT990: Bohr inger et alii 19941 : Paper X). A first effort 
towards studying the molecular gas properties of ETGs with multi- 
ple molecular t ransitions (including 13 CO, HCN and HCO + in ad- 
dition to 12 CO: ICrocker et al.l20ll hereafter Paper XI) shows that 
the molecular gas properties widely vary amoung ETGs, probably 
due to 12 CO optical thickness variations linked to the dynamical 
state of the gas. 

It is also not well determined whether the properties of in- 
dividual galaxies are the determining factors in the behavior of 
the molecular gas, or whether the molecular gas in ETGs fol- 
lows t he same star-forma ti on relations as in spi r als (e.g.lKennicuttl 
19981: iBigiel et all 120081 : IShapiro et all l20ld : IWei et all l2010al : 



Crocker et all 1201 ll) . To completely understand star formation in 



ETGs, it is vital to study in a spatially-resolved way the molecular 
gas in an unbiased sample of molecular-gas rich ETGs. Therefore, 
as part of the ATLAS 3D project, we performed the present imag- 
ing survey with the Combined Array for Research in Millimeter 
Astronomy (CARMA), the first of its kind capable of providing 
significant conclusions about the state of molecular gas in ordinary 
ETGs from interferometric imaging. 

In 33 we describe the selection criteria used to define the 
CARMA sample. In ^3] we describe the observational strategy, the 
data calibration and reduction, and the analysis techniques. In §4\ 
we investigate the properties of the molecular gas in the sample, in- 
cluding the range of morphologies present and their relationships to 
other host galaxy properties. We conclude briefly in §5\ Appendix 
lAl describes the CO data of individual galaxies in the CARMA 
ATLAS 3D sample. Appendix [B]presents CO literature data of CO- 
imaged ATLAS 3D galaxies, that will be included in the official data 
release in a uniform manner. 



2 THE CARMA SAMPLE 

The sample of ETGs chosen for the CARMA survey was extracted 
from the ATLAS survey. All ATLAS 3D galaxies were observed 
wi th the IRAM 30m te lescope in CO(l-O) and CO(2-l), mostly 
bv lCombes et al.l d2007h and Paper IV, though a few were observed 
bv lWelch & Sagel <2003l) . The CARMA sample consists of 30 of 
the CO-detected ATLAS 3D galaxies (see Table©. NGC 2697 and 
NGC 4292, removed from the ATLAS 3D sample because of evi- 
dence of spiral structure in optical images and of a lack of SAURON 
data, respectively, were also imaged with CARMA. We present 
their maps in Appendix [Cj but will not use them in any statistical 
work discussed in g] We will refer to the 30 CO-rich ATLAS 3D 
galaxies with CARMA data as the CARMA ATLAS 3D sample. 

Including 10 CO-detected sample galaxies that have already 
been interferometrically-mapped (see $\3.5\ . The total number of 
ATLAS 3D galaxies with interferometric observations discussed in 
this paper is 40. This sample of CO-imaged ATLAS 3D galaxies 
is complete for the 33 brightest objects, down to an integrated 
CO(1-0) flux limit of 18.5 Jy km s _1 (corresponding to the flux 
of NGC 3182). Two galaxies with fluxes below this threshold also 
have CARMA CO images, UGC 05408 and NGC 5173. Five other 
CO-faint galaxies h ave CO maps in the literature (NGC 4550 by 
Crocker et all | 2009 [ : NGC 2768, N GC 4477 and NGC 3489 b y 
Crocker et all l201ll : NGC 2685 bv ISchinnerer & Scovilldl2002h . 
The faintest detections from Paper IV were not observed due to 
observing time restrictions. 

Objects in the CARMA ATLAS 3D sample have total absolute 
if -band magnitudes ranging from Mr = —21.57 to —25.09, and 



distances from 13.4 Mpc to 45.8 Six of the 30 galaxies be- 

long to the Virgo cluster. All galaxies but three (IC 719, NGC 1222 
and N GC 7465) are regular ro tators according to their stellar kine- 
matics jKrainovic et al.ll201lL hereafter Paper II). NGC 1222, clas- 
sified as a non-regular rotator, is known to be undergoing a strong 
interaction with a neighbour. NGC 7465 is undergoing an inter- 
action with NGC 7464 and NGC 7463, and has a kinematically- 
decoupled core. IC 719 does show ordered rotation in its stellar 
component, but is recognized as a 2a (double-peaked) galaxy when 
the velocity dispersion map is considered, indicating two counter- 
rotating stellar discs. Further properties of the sample galaxies are 
listed in Table Q] 



3 OBSERVATIONS, CALIBRATION AND DATA 
REDUCTION 

3.1 CARMA observations 

Observations for this survey were taken in the 12 CO(1-0) line at 
CARMA over the course of five semesters, beginning in Autumn 
2008 and finishing in November 2010. Galaxies were always first 
observed with the CARMA D array, with 1 1 - 150m baselines, cor- 
responding to observable angular scales of 3.5 - 48" at CO(1-0). 
NGC 2697 and NGC 7465, that appeared to have significant flux 
resolved out, were followed up with the more compact E array (8 
- 66m baselines). NGC 1266 was spatially unresolved in D array, 
and was followed up with the CARMA B and A arrays (a de tailed 
discussion of NGC 1266 can be found in lAlatalo et al.ll20Tlh . The 
CARMA ATLAS 3D survey had a surface brightness sensitivity (la 
in 100 km s _1 line width) ranging from 12 to 369 M© pc -2 , with 
median a sensitivity of 84 M© pc -2 . On average, ^100 hours of 
observations were obtained each semester. Observational parame- 
ters are listed in Table [2] 

Upgrades to the CARMA correlator and receivers were tak- 
ing place while the ATLAS 3D survey was ongoing, so data taken 
later in the programme have larger bandwidths and simultane- 
ous observations of 12 CO(1-0) and 13 CO(1-0). A handful of 
ATLAS 3D galaxies thus also have 13 CO maps that will be pre- 
sented in an upcoming paper. The galaxies with large line widths 
(Av > 420 km s" 1 ) were observed only after the CARMA corre- 
lator was upgraded, providing sufficient bandwidth and channel 
resolution to properly cover and sample the line. We were able 
to reliably image 3mm continuum in three sources (NGC 1266, 
NGC 3665 and NGC 5866), and Table Enlists those fluxes as well 
as upper limits for the other galaxies. 



3.2 Calibration and imaging 

Raw CARMA visibility data were reduced in the usual way, using 
the Multichannel Image Reconstruction Im age Analysis and Dis- 
play (MIRIAD) package dSault et al.lll995l) . For each source and 
track, the raw data were first Hanning- smoothed in velocity. Then 
the phase vs. time behaviour of the calibrator was checked to search 
for decorrelations over baselines, which were flagged out. Next, 
the data were corrected for differences in the lengths of the fiber 
optic lines between the antennas and the correlator. The bandpass 



2 UGC 05408 was kept in the ATLAS 3D survey because its estimated dis- 
tance is within 42 Mpc when errors in the distance estimate are taken into 
account (Paper I) 
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Table 1. CARMA ATLAS 3 D galaxy sample. 



Name 


a. 


d 


^sys 




0mol 


d 


Virgo 


(J2000) 


(J2000) 


(km s _1 ) 


(cleg) 


(deg) 


(Mpc) 


membership 


IC 676 


11 12 39.84 


+09 03 20.7 


1429 


69 


16.5 


24.6 





IC719 


11 40 18.52 


+09 00 35.6 


1833 


74 


229.0 


29.4 





IC 1024 


14 31 27.07 


+03 00 30.0 


1479 


72 


24.5 


24.2 





NGC 1222 


03 08 56.76 


-02 57 19.3 


2422 


41 


33.0 C 


33.3 





NGC 1266 


03 16 00.79 


-02 25 38.6 


2170 


26 


270 


29.9 





NGC 2764 


09 08 17.44 


+21 26 35.8 


2706 


65 


202.5 


39.6 





NGC 2824 


09 19 02.22 


+26 16 12.3 


2758 


61 


161.5 


40.7 





NGC 3182 


10 19 33.02 


+58 12 21.0 


2118 


35 


331.5 


34.0 





NGC 3607 


11 16 54.54 


+ 18 03 07.1 


942 


34 


302.5 


22.2 





NGC 3619 


11 19 21.60 


+57 45 28.3 


1560 


48 


74.5 


26.8 





NGC 3626 


11 20 03.78 


+ 18 21 25.6 


1486 


67 


169.5 


19.5 





NGC 3665 


11 24 43.64 


+38 45 46.2 


2069 


64 


219.5 


33.1 





NGC 41 19 


12 08 09.60 


+ 10 22 44.7 


1656 


69 


296.0 


16.5 


1 


NGC 4150 


12 10 33.65 


+30 24 05.5 


208 


54 


146.0 


13.4 





NGC 4324 


12 23 06.17 


+05 15 02.8 


1665 


62 


232.0 


16.5 


1 


NGC 4429 


12 27 26.56 


+ 11 06 27.3 


1104 


68 


82.0 


16.5 


1 


NGC 4435 


12 27 40.49 


+ 13 04 44.3 


791 


52 


201.0 


16.7 


1 


NGC 4694 


12 48 15.10 


+ 10 59 01.3 


1171 


69 


155.5 


16.5 


1 


NGC 4710 


12 49 39.36 


+ 15 10 11.7 


1102 


86 


207.0 


16.5 


1 


NGC 4753 


12 52 22.07 


-01 11 57.9 


1163 


75 


93.0 


22.9 





NGC 5173 


13 28 25.29 


+46 35 29.9 


2424 


24 


100 c 


38.4 





NGC 5379 


13 55 34.35 


+59 44 34.3 


1774 


64 


66 


30.0 





NGC 5866 


15 06 29.60 


+55 45 48.0 


755 


89 


127 


14.9 





NGC 6014 


15 55 57.39 


+05 55 54.7 


2381 


22 


139.5 


35.8 





NGC 7465 


23 02 00.96 


+ 15 57 53.3 


1960 


70 


106.0 


29.3 





PGC 029321 


10 05 51.18 


+ 12 57 40.7 


2816 


38 


76 


40.9 





PGC 058114 


16 26 04.29 


+02 54 23.6 


1507 


71 


94.5 


23.8 





UGC 05408 


10 03 51.86 


+59 26 10.2 


2998 


31 


315 c 


45.8 





UGC 06176 


11 07 24.68 


+21 39 25.6 


2677 


68 


201.0 


40.1 





UGC 09519 


14 46 21.12 


+34 22 14.2 


1631 


41 


177.5 


27.6 





NGC 2697 a 


08 54 59.40 


-02 59 15.2 


1814 b 


30 


301.5 


22.0 





NGC 4292 a 


12 21 16.49 


+04 35 44.3 


2258 6 


50 


230 


29.8 






a Removed from the ATLAS 3D sample, but the CO data are presented in Appendix [C] 
6 From the HYPERLEDA catalog, due to the absence of SAURON data. 
c Not originally in Paper V. 
Notes: 

Column (1): Principle designation from LEDA, used as the standard designation. 

Column (2): Right ascension (J2000), from LEDA. 

Column (3): Declination (J2000), from LEDA. 

Column (4): Optical stellar heliocentric velocity, from Paper I. 

Column (5): Best-determined CO inclination angle, from Paper V. 

Column (6): Kinematic position angle of the molecular gas, from Paper X. 

Column (7): Average distance to each galaxy, from Paper I, orig inally from Tonrv et al. ( 2001) and Mei et al. ( 2007) 
Column (8): Virgo membership of each galaxy, from Paper I. 
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Table 2. Observational parameters of CARMA galaxies 



Name 


Semester 


Gain 


Total 




Ay 


Bandwidth 


calibrator 


hours 


(arcsec) 


(km s -1 ) 


(kms -1 ) 


IC 676 


2009B 


1058+015 


3.75 


3.8 x 3.3 


10 


410 


IC 719 


201 OA 


3C273 


13.17 


3.9 x 2.8 


20 


580 


IC 1024 


2008B 


3C279 


4.68 


3.9 x 3.0 


10 


420 


NGC 1222 


2008B 


0339-017 


5.20 


3.6 x 3.2 


5 


420 


NGC 1266 


2008B 


0339-017 


5.43 


4.2 x 3.3 


10 


410 


NGC 2764 


2009B 


0854+201 


4.43 


3.5 x 2.9 


10 


410 


NGC 2824 


2008B 


0956+252 


8.79 


4.3 x 4.0 


25 


400 


NGC 3182 


2009A 


0927+390 


7.42 


5.2 x 4.1 


30 


390 


NGC 3607 


201 OB 


1058+015 


9.30 


5.6 x 5.0 


20 


900 


NGC 3619 


201 0B 


0958+655 


24.44 


4.4 x 3.9 


40 


880 


NGC 3626 


2009B 


1159+292 


7.50 


3.9 x 3.7 


25 


400 


NGC 3665 


201 0B 


1159+292 


16.59 


4.3 x 4.2 


10 


910 


NGC 41 19 


2008B 


3C273 


3.67 


5.0 x 4.0 


10 


410 


NGC 4150 


201 OA 


1159+292 


6.74 


5.3 x 4.1 


10 


560 


NGC 4324 


2009A 


3C273 


20.18 


4.7 x 3.9 


20 


400 


NGC 4429 


201 0B 


3C273 


8.83 


4.7 x 3.7 


10 


840 


NGC 4435 


201 0B 


3C273 


8.48 


3.9 x 3.4 


10 


920 


NGC 4694 


201 OA 


3C273 


5.38 


3.9 x 3.1 


10 


390 


NGC 4710 


2009A 


3C273 


5.10 


3.9 x 3.2 


10 


410 


NGC 4753 


201 0B 


3C273 


10.38 


5.6 x 4.1 


15 


735 


M/^f" 1 < 1 11 

PMLrL. Jl Id 


ZUU9r> 


131U+3Z3 


1 HA 


o.y X o.O 


on 


O/i n 
Z4U 


NGC 5379 


201 OA 


1642+689 


5.33 


5.2 x 3.6 


10 


410 


NGC 5866 


201 0B 


1419+543 


3.52 


3.6 x 3.1 


10 


920 


NGC 6014 


2008B 


1751+096 


3.85 


4.2 x 3.9 


10 


350 


NGC 7465 


2008B 
2009A b 


3C454.3 


15.02 


6.6 x 5.6 


10 


410 


PGC 029321 


2009B 


1058+015 


4.52 


3.8 x 3.7 


10 


410 


PGC 058114 


2009A 


1549+026 


3.62 


4.4 x 3.8 


10 


410 


UGC 05408 


2009A 


0927+390 


2.00 


4.6 x 3.5 


25 


196 


UGC 06176 


2009A 


1159+292 


4.08 


3.5 x 2.8 


15 


735 


UGC 09519 


2008B 


1310+323 


3.75 


3.4 x 2.9 


10 


410 


NGC 2697 


2008B 
2009A b 


0927+390 


14.33 


6.8 x 5.1 


20 


380 


NGC 4292 


2008B 


3C273 


19.86 


4.1 x 3.3 


10 


410 



a Observed both in C- and D-arrays 

6 Observed in D-array in 2008B and E-array in 2009A. 

Notes: 

Column (1): Rrinciple designation from LEDA, used as the standard designation. 

Column (2): Semester in which the observations were taken at CARMA. 

Column (3): Gain calibrator used. 

Column (4): Total on-source observing time in hours. 

Column (5): FWHM of the major and minor axes of the synthesized beam. 

Column (6): Final velocity channel width. 

Column (7): Velocity bandwidth covering the line. 
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was then determined using a high signal-to-noise ratio (S/N) ob- 
servation of a bright calibrator, and used to correct for instrumen- 
tal spectral fluctuations. The atmospheric phase offsets present in 
the data were determined using a phase calibrator, usually a quasar 
within 20° of the source, observed at regular (^15 minutes) in- 
tervals. Amplitude calibration was performed using the phase cal- 
ibrator. A catalog of the most up-to-date fluxes of each calibrator 
is maintained at CARMA, both through monitoring of calibrator 
fluxes during science tracks and through dedicated flux calibration 
tracks taken on a weekly basis, and is updated weekly using all 
tracks that include a primary calib rator (e.g. a planet), used to in- 
fer the secondary calibrator fluxes dBauermeister et al]|2012h . Flux 
calibration uncertainties are assumed to be 20%, which is standard 
for millimeter observations using planetary models (Petric et al., in 
prep). After the data were satisfactorily processed, the gain solu- 
tions derived from the nearby calibrator were applied to the source. 
The calibrated source data of all observations of each source were 
then combined into one visibility file for imaging. 

We then used MIRIAD to convert all visibility files to three- 
dimensional (3D) data cubes. When possible, a zeroth order con- 
tinuum fit to the i^-data was made using channels free of line 
emission, and was subtracted from all channels. Only galaxies 
observed after the correlator was upgraded (bandwidth > 420 
km s _1 ) had sufficient bandwidth for a reliable continuum fit. 
Galaxies with narrower bandwidths were not continuum- subtracted 
(continuum emission generally adds < 10% extra flux to each 
galaxy, much less than the total flux calibration uncertainties, and 
this is usually confined to the central beam in each galaxy). The 
continuum- subtracted i^-datasets were then transformed into RA- 
Dec-velocity space (with velocities determined from the line be- 
ing imaged, mainly 12 CO but also occasionally 13 CO). Channel 
widths were chosen to achieve at least a 3 a detection in each chan- 
nel where flux was present, and were always larger (by at least a 
factor of 2) than the original spectral resolution of the i^v-data. Pix- 
els of l"x \" were chosen as a compromise between spatial sam- 
pling and resolution, typically giving approximately 4 pixels across 
the beam major axis. A constant pixel size was chosen to maintain 
uniform spatial sampling across the sample. One arcsecond corre- 
sponds to a physical scale between 72 and 222 pc, depending on 
the distance of the source. The areas imaged were generally cho- 
sen to be within the primary beam of the 10m antennas 54"), 
but in the cases of sources that extend beyond this (NGC 4324, 
NGC 4710 and NGC 5866), the area imaged was taken to be closer 
to the primary beam of the 6m antennas 90"). In these cases, the 
MIRIAD imaging task INVERT was run with the mosaicking op- 
tion, to properly scale the data and account for the different primary 
beam widths. 

Robust = weighting was used by default, but was changed 
if the size and make-up of the source dictated it. Natural weighting 
(robust = 1) was used for sources where large scale features 
presented in the initial (robust = 0) channel maps, and we sus- 
pected that some flux might be resolved out (NGC 3619, NGC 3626 
and NGC 4324). The dirty cubes were cleaned in regions of source 
emission to a threshold equal to the rms of the dirty channels. The 
clean components were then added back and re-convolved using a 
Gaussian beam of full-width at half maximum (FWHM) equal to 
that of the dirty beam. This produced the final, reduced and fully 
calibrated data cube for each galaxy. 

The 3D data cubes were then used to create moment maps of 
each galaxy: a zeroth moment (momentO) or integrated intensity 
map, and a first moment (moment 1) or mean velocity map. The 
data cubes were first Hanning- smoothed in velocity and Gaussian- 



smoothed spatially (with a FWHM equal to that of the beam), and 
masks were created by selecting all pixels above fixed flux thresh- 
olds, adjusted to recover as much flux as possible in the moment 
maps while minimising the noise (generally about 2-3 times the 
rms noise in the smoothed channels). The moment maps were then 
created using the unsmoothed cubes within the masked regions 
only. 

All data products were converted into Flexible Image Trans- 
port System (FITS) files using MIRIAD. Tables and [3] list some 
datacube properties for each of the ATLAS 3D CARMA galax- 
ies. Figure [T] shows the CO momentO contours of each CARMA 
ATLAS 3D sample galaxy overlaid on an r-band image, either from 
the Sloan Digital Sky Survey (SDSS) or the Isaac Newton Tele- 
scope (INT) when SDSS imaging is unavailable (Scott et al. 2012, 
submitted, hereafter Paper XXI). Figure [2 displays the momentl 
maps of the 30 sample galaxies. Figure [3] shows the CO momentO 
contours overlaid on unsharp-masked optical images, to enhance 
dust features (see 34. U . Figure 01 displays position- velocity dia- 
grams (PVDs) of galaxies whose kinematics are not well repre- 
sented with momentl maps, due to multiple velocity components 
along the line-of-sight (NGC 2764, NGC 4710 and NGC 5866). 
AppendixlAlcontains multiple figures for each CARMA ATLAS 3D 
galaxy, including a r-band-momentO overlay, momentO map, mo- 
mentl map, a comparison of the CARMA and IRAM 30m spectra, 
a PVD along the kinematic position angle determined in Paper X, 
and individual channel maps. A more detailed analysis of the PVD, 
as well as a comparison to the stellar and ionized gas kinematics is 
available in Davis et al. (2012, hereafter Paper XIV). 

The continuum measurement or upper limit for each galaxy 
depended on the correlator configuration. For narrow-band (band- 
width « 420 km s _1 ) galaxies that did not include 13 CO in the 
lower sideband, the full lower sideband was used to calculate the 
continuum (bandwidth ^186 MHz). For narrow-band galaxies that 
included 13 CO observations, the edges of the band in the lower 
sideband were used (bandwidth ^124 MHz). Galaxies that were 
observed with the upgraded correlator were treated individually, the 
continuum being modeled in line-free channels only, which varied 
from galaxy to galaxy. The continuum upper limit was taken to be 
three times the rms noise of the dirty map produced by inverting the 
combined line-free channels. Galaxies where a continuum source 
was detected (NGC 1266, NGC 3665 and NGC 5866), the flux was 
measured from the detected point source, and the rms noise was 
measured in the flux-free regions of the cleaned maps. 

3.3 Measurements from the data 

To calculate the noise level, cr rm s of each data cube, we selected 
all the pixels within each cube outside the region known to contain 
flux, and calculated the standard deviation. The total flux per chan- 
nel was then calculated by using the region of the momentO map 
that contained flux, then summing over all pixels in each channel 
located within that region. The statistical noise per channel is calcu- 
lated as a r ms * VNb , where Nb is the number of beam areas in the 
region with identified emission. The CARMA integrated spectra 
presented in Appendix lAlhave been constructed using this method. 
The integrated fluxes in the channel maps were converted from the 
native Jansky beam -1 units into Jansky by dividing out the total 
beam area in pixels, and are given in Table [3] 

PVDs were created by taking thin (5 pixel) slices through the 
data cubes, positioned to intersect the centre of the CO emission, 
and using the kinematic position angle determined in Paper X. Ta- 
ble Q] also lists the kinematic position angles of galaxies not orig- 
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Table 3. Data parameters of CARMA ATLAS 3D galaxies 



Name 


(Trms 

(mJy bin -1 ) 


Continuum 
(mJy bm _1 ) 


^CARMA 

(Jykms- 1 ) 


7 M Ho 


-FqARMA 
-^30m 


Morph. 
class. 


IC 676 


30.6 


< 3.09 


66.3 ±2.5 


8.71 ±0.04 


1.2 ±0.05 


D 


IC 719 


6.63 


< 0.351 


20.4 ± 1.0 


8.34 ± 0.05 


1.2 ± 0.12 


D 


IC 1024 


19.4 


< 3.06 


94.6 ± 1.7 


8.87 ± 0.02 


1.8 ± 0.04 


D 


NGC 1222 


23.1 


< 1.97 


147.6 ± 1.6 


9.33 ± 0.01 


1.8 ± 0.04 


X 


NGC 1266 


22.8 


7.49 ± 1.96 162.9 ±0.7 


9.28 ± 0.01 


1.0 ± 0.03 


D 


NGC 2764 


12.9 


< 1.85 


85.4 ± 1.5 


9.23 ± 0.02 


1.1 ± 0.04 


B+R 


NGC 2824 


10.2 


< 2.21 


37.8 ±0.8 


8.91 ± 0.02 


1.8 ± 0.08 


D 


NGC 3182 


8.16 


< 1.68 


16.4 ±0.7 


8.41 ± 0.04 


1.2 ± 0.12 


R 


NGC 3607 


7.50 


< 0.830 


58.0 ±0.8 


8.50 ± 0.01 


1.2 ± 0.04 


D 


NGC 3619 


3.23 


< 0.327 


15.2 ±0.4 


8.28 ± 0.03 


0.7 ± 0.12 


M 


NGC 3626 


5.34 


< 1.48 


43.2 ±0.9 


8.32 ± 0.02 


1.3 ± 0.10 


D 


NGC 3665 


10.2 


8.98 ±0.42 94.3 ± 1.2 


9.11 ±0.01 


1.6 ±0.06 


D 


NGC 4119 


20.2 


<2.64 


39.3 ± 0.9 


8.14 ± 0.02 


1.8 ±0.07 


R 


NGC 4150 


9.41 


< 1.87 


20.5 ± 0.4 


7.82 ± 0.02 


0.7 ± 0.08 


M 


NGC 4324 


4.30 


< 0.974 


27.2 ± 0.6 


7.97 ± 0.02 


1.9 ± 0.12 


R 


NGC 4429 


9.53 


< 0.375 


70.8 ± 0.9 


8.39 ± 0.01 


2.2 ± 0.08 


D 


NGC 4435 


11.5 


< 0.375 


31.6 ± 1.0 


8.05 ± 0.03 


1.5 ± 0.10 


R 


NGC 4694 


9.30 


< 0.512 


21.7 ±0.6 


8.01 ± 0.03 


0.7 ± 0.07 


X 


NGC 4710 


14.0 


< 5.20 


351.2 ± 1.7 


9.08 ± 0.01 


2.3 ± 0.03 


B+R 


NGC 4753 


11.1 


< 0.521 


74.7 ± 1.1 


8.70 ± 0.01 


1.4 ±0.06 


M 


NGC 5173 


8.16 


< 1.46 


12.5 ±0.2 


8.36 ± 0.02 


1.2 ±0.14 


X 


NGC 5379 


10.8 


< 0.450 


18.5 ±0.8 


8.53 ± 0.04 


1.6 ±0.15 


B+R 


NGC 5866 


17.6 


4.02 ± 0.54 258.2 ± 3.4 


8.75 ± 0.01 


1.9 ±0.02 


B+R 


NGC 6014 


22.4 


< 3.61 


34.5 ± 1.0 


8.77 ± 0.03 


1.0 ±0.07 


S 


NGC 7465 


9.71 


< 1.39 


93.6 ± 0.4 


9.02 ± 0.01 


1.7 ±0.04 


M 


PGC 029321 


13.1 


<2.09 


18.3 ±0.7 


8.61 ± 0.04 


1.2 ±0.09 


D 


PGC 058114 


15.7 


< 2.21 


73.7 ±0.9 


8.75 ± 0.01 


1.4 ±0.04 


D 


UGC 05408 


13.6 


<2.83 


6.1 ±0.5 


8.32 ± 0.08 


0.8 ±0.18 


D 


UGC 06176 


10.8 


< 0.524 


28.5 ±0.8 


8.76 ± 0.03 


1.5 ±0.10 


D 


UGC 09519 


19.5 


< 2.27 


51.4 ± 1.4 


8.77 ± 0.03 


0.9 ± 0.04 


M 


NGC 2697 


9.89 


< 1.36 


48.2 ± 0.7 


8.61 ± 0.01 


2.3 ±0.11 


R 


NGC 4292 


9.58 


< 1.05 


13.8 ±0.5 


7.74 ± 0.04 


1.2 ±0.14 


R 



Notes: 

Column (1): Principle designation from LEDA, used as the standard designation. 

Column (2): Root mean square noise level per channel, measured in regions devoid of line emission. 

Column (3): Continuum emission or upper limit (see text). 

Column (4): Total CO flux from CARMA. 

Column (5): Total molecular gas mass from the CARMA observations using Jco = 3 x 10 20 (K km s -1 ) -1 cm -2 (as in Paper IV) 
Column (6): Fraction of CO flux recovered by CARMA (with respect to the IRAM 30m telescope). 
Column (7): CO morphological classification, based on the momentO and momentl maps: 
D = disc, X = strongly disrupted, R = ring, M = mildly disrupted, B = bar/ring, S=spiral. 



inally listed in Paper X, but subsequently had kinematic position 
angles fitted using the same methods. 

3.4 Comparison to the 30m data 

The integrated spectra from CARMA and the IRAM 30m telescope 
are overlaid in Figs. IA1I4A30I Millimeter measurements are as- 
sumed to have 20% absolute flux uncertainties (errors reported on 
the measured parameters listed in Table [3] only include the random 
errors on the measurement itself, not the systematic flux calibration 
error), so flux ratios between 0.78 and 1.28 indicate fluxes in agree- 
ment within the uncertainties. The spectra from the IRAM 30m 
agree with the CARMA ones for 15/30 galaxies. In all cases where 
there is disagreement, CARMA actually recovers more flux than 
the 30m. Most discrepant galaxies either extend beyond, or fill the 
majority of the single-dish beam, with two large ratios (NGC 4435 
and UGC 06176) likely due to a slight off-center pointing of the 
single dish. In fact, only 7 of the galaxies within the CARMA 
ATLAS 3D sample are well contained within the 21.6" IRAM beam. 



As discussed above, all galaxies detected with the IRAM 30m 
were also detected at CARMA. However, the IRAM 30m selec- 
tion for the current CARMA survey does imply that galaxies with 
molecular gas exclusively beyond the 30m beam (e.g. rings, arms 
and tidal features beyond a radius of w 12") would have been 
missed. It is impossible to estimate the incidence of these objects 
based on the CO data alone, but the Galaxy Evolution Explorer 
(GALEX) ultraviolet imaging survey of the SAURON sample (trac- 
ing the associated star formation) suggests that they a re very rare 
(of 34 ETGs, only one has such features, NGC 2974; Ijeong et all 
l2Q09h . Lower density neutral hydrogen is however fairly common 
at large radii (Paper XIII). 



Because the molecular gas in many (> 75%) galaxies extends 
beyond the single-dish beam, it is likely that the CARMA flux more 
accurately reflects the total molecular gas mass in each galaxy, and 
therefore we adopt the CARMA derived molecular gas masses for 
the remainder of this paper. 
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Figure 1. r-band images (grayscale and black contours) overlaid with integrated CO(1-0) contours (yellow) at [10%, 30%, 50%, 70%, 90%] of the peak 
intensity for the 30 galaxies in the ATLAS 3D CARMA sample. The r-band images were either taken from the SDSS or a dedicated program with the INT 
(Paper XXI) The CO synthesized beam is shown in the bottom-right corner of each panel. A bar in the bottom left corner of each panel indicates the scale of 
10", as well as the equivalent physical scale at the distance of the galaxy. 
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Figure 2. Mean CO(1-0) velocity maps of the 30 galaxies of the ATLAS 3D CARMA sample. Iso-velocity contours (black) are overlaid at 20 km s _1 intervals, 
as are contours from the momentO map (white). The CO synthesized beam is shown at the bottom-right corner of each panel. A bar in the bottom left corner 
of each panel indicates the scale of 5 /x , as well as the equivalent physical scale at the distance of the galaxy. 
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3.5 Literature data 

Twelve galaxies of the ATLAS 3D sample already have in- 
terferometric CO data av ailable in the lite rature. Those 
galaxies are: NGC 4 710 dWrobel et al.l Il992h . NGC 2685 
dSchinnerer & Scovilld l2002h . and the ten galaxies th at will 
be p art of the ATLAS 3D data release: NGC 4476 (lYound 
l2002h: NGC 3032, NGC 4150, NGC 4459 an d NGC 4526 
(lYoung. Bureau & Cappeliaril l2008n: N GC 2768 dCrocker et all 
l2008h : NGC 455 dCrockeret alJl2009l) : NGC 524, NGC 3489 
and NGC 4477 dCrocker et all 1201 lb . Names, distances, Virgo 
membership, CO morphological classes and references of these 
galaxies are listed in Table |4] When required, data for each of 
the galaxies were provided by the original authors for use in this 
work. Moment maps and dust comparisons appear in Figure IB 1 1 
for the galaxies whose data will be included in the data release. 
Both NGC 4710 and NGC 4150 were re-observed as part of the 
CARMA ATLAS 3D survey, and we only present these updated 
data in this paper. 



3.6 Catalog of ATLAS 3D interferometric CO data 

The CARMA data archival will include CARMA data for all 30 
galaxies, as well as data for the 10 previously published galaxies 
and the two non- ATLAS 3D sources. The data products provided 
will include data cubes and momentO and moment 1 maps in FITS 
format. The data will also be reachable from the ATLAS 3D project 
websitjf). 



4 DISCUSSION 

Figures Q] and [2 show the 30 ETGs imaged as part of the 
CARMA ATLAS 3D survey. Before the survey was complete 
(which includes the 10 literature galaxies mentioned in ^3.5K 



Table 4. Properties of the ATLAS 3D literature galaxies 



16 ETGs had spatially-resolved CO maps (jWrobel et al 
Schinnerer & Scovilld I2OO2I : lYoung 2002, 2005; W ei et al.l 



[1992 



2010a; 



Young, Bureau & Cappellari 2008; Crocker et al. 2011), most ob 



tained as part of individual or small- sample studies. The addi- 
tion of the CO-imaged ATLAS 3D sources increases this amount 
to over 50, and provides a largely unbiased sample. Combined with 
the maps from the SAURON integral-field spectrograph, the CO- 
mapped ATLAS 3D galaxies thus provide the most robust picture of 
molecular gas-rich ETGs. 

We now focus our analysis on the CO-rich ATLAS 3D ETGs 
only, adding the ten ATLAS 3D galaxies imaged in CO with the 
Berkeley-Illinois-Maryland Array (BIMA), the Owens Valley Ra- 
dio Observatory (OVRO) and the Plateau de Bure Interferometer 
(PdBI) prior to the CARMA effort. Between the 30 CARMA galax- 
ies and 10 previously published galaxies (excluding two overlap- 
ping objects; see Table |4]), we imaged 40 galaxies or 65% of the 
CO detections in Paper IV, more than 90% of the total molecular 
gas mass in those detections, and are complete down to a CO(1-0) 
integrated flux of 18.5 Jy km s _1 . 



4.1 Comparison to the dust 

Figure [3] shows the CO momentO contours of each galaxy overlaid 
on — r colour image. The dust distributions are shown by creat- 



3 http://carma-server.ncsa.uiuc.edu:8181/ 

4 http://purl.org/atlas3d 



Name 



Distance Morph. 



NGC 524 
NGC 2685 
NGC 2768 
NGC 3032 
NGC 3489 
NGC 4459 
NGC 4476 
NGC 4477 
NGC 4526 
NGC 4550 



(Mpc) 


class. 


23.3 


D 


16.7 


R 


21.8 


D 


21.4 


D 


11.7 


S 


16.1 


D 


17.6 


D 


16.5 


R 


16.4 


D 


15.5 


D 



Virgo 
membership 



Reference 



Re ferences: (1) "Schinnerer & Scoville (2002); (2) Young (2002); 
(3) lYoung. Bureau & Catroellaril d2008h : (4) [Crocker et all d2QQ8h : (5) 
ICrocker et al.N2009h : f6) ICrocker et aflfcOlll 
Notes: 

a CO morphological classification based on the momentO and momentl 
maps. D = disc, X = strongly disrupted, R = ring, M = mildly disrupted, 
B = bar/ring, S=spiral. 

ing g - r colour SDSS images (and for NGC 1222 and NGC 7465, 
for which SDSS images are unavailable, on INT images, and for 
NGC 1266, on CTIO b and r band images from the SINGS survey; 
Kennicutt et al .120031) . In order to create the g — r images, the God- 
dard IDL library task SKY was used to determine the background 
sky level, then subtracted from the images. The g — r colour im- 
ages were then created by dividing the sky- subtracted g-band image 
from the sky- subtracted r-band image. The g — r imaging enhances 
the visibility of dust features by directly tracing extinction and pro- 
vides a qualitative view of the distribution of the filamentary dust 
in each galaxy of the ATLAS 3D CARMA sample. This dust, for 
the most part, is a reasonable indicator of the distribution of the 
molecular gas. 

For simplicity we refer to the g — r colour images as "dust 
images" for the remaining text, even if they are intrinsically only 
representations of stronger dust-obscured regions. Indeed, the CO 
is nearly always (28/29 instances) associated with evidence of 
dust obscuration. It is evident that the converse is not true, as a 
Hubble Space Telescope (HST) survey of nearby elli pticals re- 
vealed that they contain dust in nearly half the cases (iTran et al] 
2001: lLaueret all 12005). much more than the CO detection rate 
of ATLAS 3D . Additional dust is seen beyond the limits of the de- 
tected CO emission in most galaxies, but notably in NGC 1266, 
NGC 3619, NGC 4150, NGC 4694 and NGC 4753. However, there 
is no noticeable dust beyond the detected extent of the CO in oth- 
ers. It thus appears that in most cases where the dust and CO are 
co- spatial, dust provides a good, high spatial resolution rendering 
of the molecular gas structure. Where dust is more extended than 
the CO, it could simply be that the outer CO is below our detec- 
tion threshold, or that the ISM is not molecular at all (for instance, 
it could be atomic or warm). With this caveat in mind, we use the 
dust images as a useful guide to the morphological structures in 
which the cold ISM resides. We exploit this fact below when the 
CO morphological class is difficult to determine (e.g. IC 676 and 
NGC 4435). 



4.2 Morphologies of the molecular gas 

The CO-imaged ATLAS 3D survey illustrates that the molecular 
gas in ETGs comes in a large variety of morphologies. We have 
identified six different morphologies, including smooth discs (D), 
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Figure 3. Comparison of the integrated CO(l-O) distribution (magenta contours) and the dust distribution as represented by g — r colour images (grayscale) 
for the 30 galaxies of the ATLAS 3D CARMA sample. It appears that galaxies with disrupted CO morphologies also exhibit disturbance in the dust distribution, 
and overall the molecular gas in the galaxies traces the filamentary dust. 
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Table 5. CO morphologies of the ATLAS 3D interferometric 
sample. 



Morphological 
configuration 


Symbol 


Number 
(out of 40) 


Fraction 

(%) 


Disc 


D 


20 


50 


Spiral 


S 


2 


5 


Ring 


R 


6 


15 


Bar+Ring 


B+R 


4 


10 


Mild disruption 


M 


5 


12.5 


Strong disruption 


X 


3 


7.5 



spirals (S), rings (R), bars+rings (B+R), as bars are always found 
with rings, never alone, mildly disrupted objects (M), and strongly 
disrupted objects (X). These classifications are strictly for the 
molecular gas, as all galaxies are optically-classified as ETGs. The 
CO morphological classification of each galaxy was based on an 
analysis of its momentO and dust maps (Fig. [3]), momentl map 
(Fig. [3 and major-axis PVD (Figs.lAll-30). In many cases the mo- 
mentO and momentl map allowed us to unambiguously classify 
a galaxy. For ambiguous cases, the PVD was used to search for 
distinct kinematic signatures. Finally, the dust map was used as a 
higher spatial resolution proxy for the molecular gas. 

The criteria for each CO morphological class are as follows. 
Galaxies classified as discs (D) show regular rotation and a smooth 
elliptical shape. The PVD of D-class objects also show smooth ve- 
locity variations with no discontinuity. Galaxies classified as spi- 
rals (S) show regular rotation, but the momentO and/or dust maps 
exhibit discrete spiral arms. Galaxies classified as rings (R), like 
discs, exhibit regular rotation, but unlike discs, the PVD of R-class 
objects include a solid-body component and/or a discontinuity in 
velocity. Well-resolved rings show a central hole in the momentO 
map. In the most ambiguous cases, the dust map was consulted to 
search for a dust ring (e.g. NGC 4119 and NGC 4435). Bar+Ring 
systems (B+R) were identified either directly from the momentl 
map, when the bar and ring are both distinctly visible, or using the 
PVD and searching for an X- shaped pattern (in the case of edge-on 
systems; see ^4.2.21 ). Galaxies classified as mildly disrupted (M) 
have mainly regular disc morphologies, but include small devia- 
tions in their momentO and momentl maps. Finally, a galaxy is 
classified as strongly disrupted (X) for showing strong irregulari- 
ties in the CO kinematics. . 



4.2. 1 Strongly and mildly disrupted objects 

The origin of the strongly disrupted (X) objects is the most easily 
explained. These galaxies show highly irregular CO distributions 
and kinematics, that are normally characteristic of interacting sys- 
tems, a fact often confirmed by observations at other wavelengths. 
In the case of NGC 1222 and NGC 4694, it is clear that they ac- 
quired their ga s via external accretion (see iBeck et aD 120071 and 
iDuc et al ]|2007l for details on each galaxy). NGC 5173 is known to 
host a large misaligned H I disk (Paper XIII). 

Mildly disrupted CO distributions (M) also appear to be 
in less relaxed systems than the rest. In some cases, there is 
strong evidence t hat the gas is being acquired externally, such 
as in NGC 7465 jLi & Seaquistlll994rh which is accreting from 
an external H I reservoir via a close interaction with other galax- 
ies. NGC 3619 and UGC 09519 are known to host large mis- 
aligned H I discs (Paper XIII), providing a likely source for the 
molecular gas. The irregular dust filaments in NGC 4753 are 
the main reason that the CO is classified as an M morphology. 



The NGC 4753 dust structure is well explained by a precessing 
twisted disc after the accretio n of a gas-rich dwarf companion 
JSteiman-Cameron et al. 1 1992). The original BIMA observations 
of NGC 4150 (lYoung. Bureau & Cappellaril l2008h show a pro- 
nounced low- surface density elongated tail that was not recovered 
by CARMA, both due to a smaller field of view and higher spa- 
tial resolution. A study of the gas-phase metallicity of NGC 4150 
also show it to be sub-Solar (Z « 0.3 — 0.5 Z©), therefore the 
molecular gas is likel y to originate from an external accretion event 
dCrockettetal.ll201ll) . 

It appears that field galaxies are far more likely to exhibit 
disrupted morphologies than Virgo cluster galaxies, with 24 ± 8% 
(7/29) of field galaxies classified as either M or X as opposed to 
9 =b 9% (1/11) for Virgo galaxies. This is not surprising given the 
results of Paper X, which suggested that galaxies virialised within 
the cluster are unable to replenish their stores of molecular gas via 
external accretion, and thus only longer-lived phenomena survive, 
such as bars and rings. Overall, the M- and X-class objects seem to 
host gas that has an external origin, and they are probably dynami- 
cally younger than the other objects, explaining the dearth of such 
objects within Virgo. 



4.2.2 Molecular rings and bars 

Rings (R) and bars+rings (B+R) also make up a non-negligible 
fraction of the ATLAS 3D molecular gas morphologies, represent- 
ing about a fifth (8/40) of the CO sample. It is of note that cold 
gas in barred galaxie s accumulates in rings associated with the ma- 
jor resonances (e.g. iButa & Combesl fl996k so it is possible that 
bar+ring systems are simply in an early stage of an evolution- 
ary sequence, whereby the bar funnels gas into the centre (e.g. 
iBournaud & Combesl l2002h . evacuating the central (but not nu- 
clear) regions of molecular gas, and later leaving the ring struc- 
ture intact, without a molecular gas-rich bar. This scenario would 
explain why bars are only seen with rings, but not the converse. 

Figure 0] shows the PVDs of three galaxies (NGC 2764, 
NGC 4710 and NGC 5866) that possess two velocity components 
along the line-of- sight. All three galaxies are edge-on, and the char- 
acteristic X-shaped PVDs observed are normally associated with 
the kinematic signature of a bar in an edge-on disc. The mate- 
rial in the inner rapidly-rising component of each PVD is asso- 
ciated with gas on X2 orbits (elongated perpendicular to the bar 
and located within the inner Lindblad resonance, ILR), while the 
material in the outer slowly-rising component is associated with 
gas on nearly circular orbits beyond corotation. The intermediate 
region has been largely s wept free of gas by the bar, due to in- 
flows toward the ILR (see lSellwood et al.ll 19931 for a review of bar 
dynamics). The signature is analogous to that of the Milky Way 
bar in a Galactic longit ude-velocity diagram. It was first proposed 
in external galaxies bvlKuiiken & Merrifieldl dl995[). later refined 
by IBureau & Athanassoulal dl999h and lAthanassoula & Bureau! 
( 1999 ), who both identified multiple examples obse rvationally 
(iMerrifield & Kuiikenl Il999l : IBureau & Freeman! Il999h . The bar 
signatures in NGC 2764 and NGC 5866 are observed here for the 
fir st time, but that in N GC 4710 was previously observed in CO 
bv lWrobel et all (fl992h. The link between bars and boxy- shaped 
bulges (such as that in NGC 4710) has long been established, but 
the presence of a bar in NGC 5866 is rather unexpected given its 
large classical bulge fsee lAthanassoulall2005l for the different bulge 
types and the bar-boxy bulge relationship). 

The fraction of molecular gas-rich ETGs (6/40) that contain a 
ring structure is also a factor of « 2 larger than that of LTGs, based 
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Figure 4. Position-velocity diagrams of the three CARMA ATLAS 3D ETGs with clear two-component velocity profiles, that are not properly represented by 
the moment 1 maps. The galaxies are NGC 2764 (left), NGC 4710 (middle) and NGC 5866 (right). All 3 galaxies are edge-on or almost edge-on and show the 
characteristic X-shaped PVD signature of an edge-on bar (see ^4.2.21 



on the BIMA Su rvey of Nearby G alaxies (SONG) sample (3/39; 
ISheth et al . 2002; H eifer et al.l2 003) modified to match the average 
noise and redshift of the ATLAS 3D sample (as in Paper XIV and 
£14.2.31 ). Also, unlike for LTGs, the molecular gas morphology does 
not necessarily follow the stellar morphology, most notably stellar 
bars. Galaxies with stellar bars identified in Paper II are as likely to 
have a disc-like (D) CO morphology as R and B+R morphologies. 

It also appears that Virgo cluster galaxies are slightly more 
likely to have R and B+R structures than field galaxies, with 
45 ± 15% (5/11) of Virgo galaxies classified as either R or B+R, as 
opposed to 17 ± 7% (5/29) for field galaxies. Paper X has shown 
that the molecular gas in ETGs is preferentially kinematically- 
aligned with the stars in Virgo galaxies, so it is perhaps not surpris- 
ing to find relatively more bars and rings there, where the gas has 
been under the influence of secular processes without external gas 
interference for some time. This may also point to these morpholo- 
gies being long-lived, as opposed to the M and X morphologies 
discussed above. It is also quite possible that the slightly higher de- 
tection rate of R and B+R structures in Virgo with respect to the 
field galaxies is due to a spatial resolution effect, as Virgo g alax- 
ies are closer on average (D « 16.5 Mpc: lMei etaIll2007h than 
the field galaxies in the sample (D « 28 Mpc). If we redshift the 
Virgo galaxies to the average distance of field galaxies, the detec- 
tion rate of R and B+R structures drops from 45 to 18% and is then 
statistically consistent with the field rate (NGC 4119, NGC 4435 
and NGC 4477 become unresolved). Higher resolution and deeper 
imaging is required to determine whether the higher rate in Virgo 
is truly a resolution effect or an intrinsic property. 



4.2.3 Molecular discs and spirals 

Smooth molecular discs make up the largest fraction of the CO 
morphologies that are seen in the CO-imaged ATLAS 3D galaxies, 
while spirals make up the smallest fraction. This relative dearth 
of discrete spirals as compared to smooth discs seems to indicate 
that there is something fundamentally different about the way in 
which molecular gas behaves within an ETG as opposed to a LTG. 
A caveat to this scenario is that galaxies with dominant molecular 
gas spirals (and associated star formation) would likely be classi- 
fied as LTGs, and thus be excluded from the sample. 

It is possible that molecular gas in ETGs remains in a smooth 
distribution and lacks distinct spiral arms for dynamical reasons. 
For example, the steep potential wells in the central parts of ETGs 
(where the molecular gas resides) yield high epicyclic frequencies, 



and may thus prevent gravitational instabiliti es linked to sp iral arms 
from developing (Toomre Q parameter; e.g. lToomre1ll98lh . 

It is also possible that some of the D morphologies would 
have been classified otherwise with observations of better sensi- 
tivity and/or spatial resolution. For instance, much higher spatial 
resolu tion CO observations of NGC 1266 presented in lAlatalo et al.l 
(201 l|) reveal a disrupted nuclear disc and a molecular outflow, both 
of which are unresolved here. In addition, if only the peak of a spi- 
ral structure were detected, and convolved with a large beam, it 
could well be c lassified as a disc. In fact, as much as half of BIMA 
SONG spirals faelfer et alJl2003h would be classified as D if they 
were located at the average distance of the ATLAS 3D sample, with 
similar noise properties and beam sizes (see Paper XIV for a full 
discussion of properly redshifting the BIMA SONG data to match 
ATLAS 3D ). However, even assuming that half of the current CO- 
rich ETGs classified as D actually harbour spatially-unresolved spi- 
rals, the spiral fraction of ATLAS 3D galaxies would be 30 ± 7%, 
still much lower than the spiral fraction of BIMA SONG galaxies 
(62 ± 10%). 

While obtaining higher resolution molecular gas maps is the 
only rigorous way to directly break the S/D degeneracy, we use here 
the unsharp-masked dust images as proxies. As mentioned in £14.11 
the molecular gas morphology follows the dust morphology faith- 
fully in the majority of cases. The unsharp-masked images there- 
fore may be utilized to provide the higher spatial resolution and 
sensitivity required to separate the S and D classifications. The dust 
comparison indicates that in the majority of cases (with the possible 
exceptions of NGC 3619, NGC 3626 and NGC 3665), and contrary 
to LTGs, a disc morphology without distinct spiral arms is indeed 
the most common morphology for the molecular gas in ETGs. 



4.2.4 Summary of the gas morphologies 

Table[5]lists the total number of ATLAS 3D ETGs within each CO 
morphological class. Overall, it thus appears that the CO in ETGs 
is most likely to be in a settled morphology, either in regular discs 
(comprising 50% of the galaxies) or ring, bar+ring and spiral mor- 
phologies (comprising another 30% of the sample). This means that 
the majority (80%) of ETGs have a dynamically settled molecu- 
lar gas morphology. The dynamically unsettled (M and X) systems 
comprise the remaining 20%, and just over half of those are only 
mildly disrupted. Of course, the CO structures observed here are 
rather compact spatially, with an average extent of only 1 kpc (Pa- 
per XIV), and have correspondingly small dynamical timescales 
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Figure 5. Color-magnitude diagram of the CO-imaged ATLAS 3D galaxies, 
sorted based on the CO morphological clas sifications. The con tours (black) 
are from the SDSS Data Release 8 (DR8: lAihara et alJl201 ll) . and repre- 
sent all galaxies with redshift z > 0.08. Data for the ATLAS 3 D ETGs 
are either from the SDSS DR8 u, r and M r catalog products or the INT 
(NGC 524, NGC 3489 and NGC 7465; Paper XXI). Colour information is 
unavailable for NGC 1222 and NGC 1266 due to poor u-band photometry, 
for PGC 0581 14 because of the presence of a bright star less than a degree 
from the galaxy, and for NGC 4710 due to erroneous SDSS photometry. 



(10 7 — 10 8 yrs). The CO distributions are thus expected to relax 
and reach equilibrium rapidly. 

Figure \5\ shows the u — r vs. M r color-magnitude diagram 
(CMD) of the CO-imaged ATLAS 3D sample, from SDSS and INT 
data, overlaid upon the SDSS D ata Release 8 distr ibution of all 
galaxies with redshift z > 0.08 (lAihara et alJl201lh . There is no 
clear separation of the morphological groups, although some trends 
do appear. The brightest objects (M r < — 21) tend to be classified 
as D galaxies, with the exception of NGC 4753, which was dis- 
cussed in £14.2.11 The CO discs in these massive galaxies are also 
generally well-resolved, limiting the chance that these objects are 
poorly resolved rings or spirals. It also appears that although B+R 
galaxies can be found throughout the CMD, galaxies classified as 
R only seem to cluster on the red sequence, perhaps indicating a 
later evolutionary stage. 



4.3 CO spectral asymmetries 

It is worth mentioning that a fraction of the CO spectra exhibit 
asymmetries, and have representatives in many CO morphologi- 
cal classes (as opposed to solely M and X-class objects). Although 
asymmetries detected from a single-dish spectrum could be due to 
a slightly offset pointing, interferometric identifications of spec- 
tral asymmetries are more robu st. Asymmetric spec tra have long 
been known within H I discs fe aid win et alj Il980h . Within the 
Westerbork H I sample of Spiral and Irregular Galaxies (WHISP) 
at least 80% of all galax ies contain some amount of asymmetry 
(Ivan Evmeren et al.ll201ll) . In fact, asymmetries are found at small 
radii in at least 60% of H I discs from the WHISP survey. Despite 
the different dynamical timescales at the radii where the CO has 
been found, finding a few asymmetric objects within the CARMA 
ATLAS 3D sample is unsurprising. 



Asymmetric H I spectra generally trace lopsided spatial distri- 
butions, and lopsided molecular gas discs are likely a result of the 
molecular gas responding to a change in the gravitational poten- 
tial. In particular, long-lived lopsided H I discs are well reproduced 
in models where the gas is off-ce ntre with respect to the galactic 
halo (e.g. iLevine & Sparkel 1 1 99 8l) . Other possible causes for lop- 
sidedness include responding to an interaction (e.g. I Jog & Combes] 
l2009h . or possibly an internal instability that perturbs both the 
m — 1 (lopsided) and m — 2 (barred) modes (iMasset & Tagger] 



ll997l : lBournaud et a l. 2005; J og & Comb es 2009). In the strongest 
case of asymmetry (IC 1024), it is likely that the molecular gas 
disc is responding to a tidal interaction, which is also seen in the 
optical data (Paper II). IC 1024, despite being lopsided, is dynami- 
cally relaxed (Paper XIV) and shows regular kinematics and a CO 
disc morphology in the momentO map. This may suggest that the 
tidal interaction happened sufficiently long ago for any detectable 
morphological and kinematic disturbance of the molecular gas to 
subside, while preserving the spectral asymmetry. 

In the context of H I, the structures are well resolved, ex- 
tended, and clearly "lopsided", and the Hi cloud size is small 
compared to the global structure, and thus inferences can be made 
in most cases that the driver of the asymmetries are global, on a 
long timescale. Structures traced by the CO distributions are much 
less extended compared to the rest of the galaxy (Paper XIV), and 
thus the interpretation of CO asymmetries is less straightforward. A 
more comprehensive study of asymmetries and lopsidedness in CO 
distributions within galaxies would likely address these issues, as 
well as uncover many other yet unstudied properties of molecular 
gas distributions in galaxies, but is beyond the scope of this paper. 



4.4 Mass of externally-acquired molecular gas in ETGs 

Paper X discusses the origins of the molecular gas in ETGs, and 
identifies galaxies that likely have externally- acquired gas. Using 
the criteria set forth in Paper X, that the kinematic major axis of 
the ionized gas must be misaligned with the kinematic major axis 
of the stars by at least 30°, at least thirteen CO-imaged ATLAS 3D 
galaxies are identified as having externally- acquired gas: IC 719, 
NGC 1266, NGC 2685, NGC 2768, NGC 3032, NGC 3626, 
NGC 4694, NGC 5173, NGC 7465, PGC 029321, PGC 058114 and 
UGC 09519. NGC 3619 is classified as M in this work, exhibiting 
an unsettled CO morphology, but the origin of the CO is difficult to 
determine, with a kinematic misalignment angle of only 24 ± 3°, 
right at the boundary of being classified as external (Paper X), 
and the galaxy contains a kinematically-misaligned H I disc (Paper 
XIII). Although it is possible that the CO indeed has an external 
origin, there is no conclusive measurement to classify it as such, 
and NGC 3619 will not be considered as having external CO in the 
following discussion. The CO of NGC 4150 is included here as be- 
ing external in origin, as studies show that its gas-phase metallicity 
is sub-Solar, and thus is very li kely from the accretio n of a gas-rich 
companion in a minor merger dCrockett et alj|201l|). NGC 4 753 is 
classified as M in this work, and ISteiman-Cameron et al] (l992) 
make a compelling argument that the molecular disc is unsettled, 
and must come from an external source. For this reason, we in- 
clude here NGC 4753 in the category of gas acquired externally. 
The disrupted stellar and gas kinematics of NGC 1222 (CO class 
X), which is curr ently undergoing a strong three-way interaction 
(Beck et al. 2003), also suggest that the molecular gas is likely to 
be of external origin. Including these three additions, 15 (38%) of 
the 40 gas-rich CARMA ETGs contain external gas, in good agree- 
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ment with the fraction of externally- acquired gas estimated in Paper 
X (36%). 

The fact that the majority of galaxies with unsettled gas mor- 
phologies are identified as having an external gas origin is un- 
surprising, as a parcel of gas falling into a galaxy takes time to 
virialise, which is not the case for internal gas (stellar mass loss). 
The rest of the galaxies with external gas have disc morphologies 
(IC 719, NGC 3032, NGC 3626, PGC 029321 and PGC 058114), 
represented at about the same rate as in the sample as a whole. 
NGC 2685 and NGC 2768 are both known polar rings/d iscs, 
and are discussed i n detail in ISchinnerer & Scovillel J2002h and 
ICrocker et all (12008). respectively. 

We use the Fb masses from CARMA, where available, and 
find that the sum of the molecular gas masses in the 15 sources with 
external CO is 8.26 x 10 9 M fi , assuming X C o = 3 x 10 20 (K km 
s -1 ) -1 cm -2 (jPickman. Snell & Schloerbl 19861) . and accounts for 
44% of the total molecular gas mass in our sample of 40 CO-rich 
galaxies with interferometric maps. There is reasonable agreement 
between the external number (38%) and mass (44%) fractions. We 
note that these fractions are lower limits, due to the fact that we 
can only unambiguously assign galaxies with misaligned stellar 
and gas kinematics as external accretors, while some fraction of 
aligned galaxies are also likely to be so. 

Galaxies that were originally detected in molecular gas (Paper 
IV) but were not mapped interferometrically have a higher rate of 
ionised gas kinematically- aligned with the stars, compared to those 
actually mapped here (i.e. brighter CO; see Paper X). This likely 
means that the un-imaged CO-detected ETGs contribute negligibly 
(< 6%) to the total external molecular gas mass. It is also likely 
that undetected galaxies do not contribute significantly to the total 
gas mass, with an upper limit to the contribution of^3xl0 9 M©, 
assuming that 42.5% upper limits in undetected field galaxies are 
due to an external phenomenon (Paper X). 

Assuming the molecular gas is depleted through star for- 
mation and reasonab le gas depletion timescales 1 Gyr; 
ISaintonge et al.l Eoill) and star formation histories (allowing us 
to adopt the mean curr ent-to-initial molec ular gas mass correc- 
tion factor calculated by iKavirai et al .1120121 for ETGs with promi- 
nent dust lanes), we can infer the amount of molecular gas each 
galaxy acquired before star formation activity consumed some of 
the molecular gas. The total amount of molecular gas originally 
acquired externally within the ATLAS 3D sample is then approxi- 
mately 2.48 x 10 10 M©. At the rate of 3 minor mergers in the 
sample detailed in £Q]and from the rate of lLotz etaD J2008I) . this 
would mean approximately 8.3 x 10 8 M© of H2 per minor m erger 
per molecular gas depletion time (1 Gvn lSaintonge et al.l201ll) . We 
derive an average stellar mass ratio for these minor mergers by first 
converting the molecular gas mass to an atomic mass using th e 1:3 
H 2 -to-Hl ratio seen in the Milky Way dWolfire et al.ll2003l) . We 
then use the aver age H I- to- st e llar m ass fraction of minor merger 
companions from iKannappanl J2004I) of 20% to calculate the av- 
erage stellar mass of the accreted companions, 1.25 x 10 10 M©. 
For the mass of the accreting galaxy, we take the average stellar 
mass of the non- Virgo ETGs in the CO-imaged ATLAS 3D sample, 
« 6.5 x 10 10 M©. This results in a merger stellar mass ratio of 
1:5.2, consistent with the minor merger stellar mass ratios we ex- 
pect. We therefore conclude that our results are in good agreement 
with the origin of the external gas being due to minor mergers. 



5 CONCLUSIONS 

We have presented millimeter-wave data products as part of the 
CARMA ATLAS 3D survey of galaxies, an imaging survey of the 
CO J = 1-0 emission of 30 nearby CO-rich ETGs, part of the 
ATLAS 3D survey. The main conclusions drawn from these data are: 

1. The CO-imaged ATLAS 3D survey is the largest CO imag- 
ing survey of ETGs to date, and provides the most detailed view of 
the nature of molecular gas in ETGs. 

2. The molecular gas found in ETGs appears to be co- spatial 
with dust features (as traced by g — r colour images), though the 
dust features often appear to be more extended, likely because we 
are unable to trace small molecular gas surface densities. 

3. The molecular gas present in ETGs comes in a variety of 
morphologies: 50% of the objects have their CO in inclined gas 
discs or is poorly resolved structures (D), 5% show spiral arms 
or spiral structure (S), 15% have resolved rings (R), 10% have 
bar+ring systems (B+R), 12.5% have mildly disrupted distributions 
(M) and 7.5% have strongly disrupted distributions (X). 

4. The CO morphology does not show a strong correlation 
with the u — r color or total magnitude of the galaxy, although there 
appear to be weak trends. For example, the most massive galaxies 
in this sample tend to have D CO morphologies, and R morpholo- 
gies tend to lie on the red sequence whereas B+R morphologies 
tend to also be present in the blue cloud. 

5. We currently observe a total of 8.26 x 10 9 M© of molecular 
gas that is likely of external ori gin in 15/40 of our g alaxies. Using 
the correction factor derived in K avirai et all (120121) . we estimate 
that the total externally acquired molecular gas prior to the con- 
sumption of some of it by star formation was « 2.48 x 10 10 M©. 
In conjunction with the minor merger rate of lLotz et"aD j2008l) and 
assumptions about gas mass fractions, this amount of externally ac- 
creted molecular gas is consistent with a minor merger origin. 

6. The CARMA ATLAS 3D (as well as 9 additional literature 
galaxy) data products including data cubes and momentO and mo- 
ment 1 maps in FITS format, will be publicly released and hosted 
through the CARMA data archival and will be reachable from our 
project websit^E 
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APPENDIX A: CO CARMA DATA OF ATLAS 3 GALAXIES 

For each galaxy of the CARMA ATLAS 3D sample, we provide six figures: 

Top-left: r-band image (grayscale) either from SDSS or from the INT, overlaid with black contours. The CO integrated intensity (momentO) 
map is also overlaid (cyan contours corresponding to 20, 40, 60, 80 and 100% of the peak). The cyan box corresponds to the size of the 
individual moment maps in the middle-left and middle-centre panels. 

Top-right: Comparison of the single-dish IRAM 30m CO(1-0) integrated flux (black) and the CARMA integrated flux (red). The total 
integrated CARMA flux is obtained by summing all flux above a given threshold (see £13. 3b . The grey dashed vertical lines indicate the total 
spectral range of CARMA. Both heliocentric velocities and velocities with respect to systemic (taken from Paper I) are indicated. The rms 
noise is shown as a single error bar in black for the IRAM 30m and red for CARMA. 

Middle-left: CO(1-0) integrated intensity (momentO) map (colorscale), overlaid with black and grey contours (at equal increments of 20% 
of the peak). The colour table on the right provides the integrated flux scale. The dot-dashed gray circle corresponds to the 21 / /6 primary 
beam of the IRAM 30m telescope. The CARMA synthesized beam is shown in the bottom-right corner (black ellipse). 

Middle-center: CO(1-0) mean velocity (momentl) map (colorscale and black contours at 20 km s _1 intervals unless otherwise stated). 
The colour table on the right provides the velocity scale with respect to the systemic velocity (taken from Paper I and indicated in the 
top-left corner of the panel). The dot-dashed gray circle corresponds to the 21 / /6 primary beam of the IRAM 30m telescope. The CARMA 
synthesized beam is shown in the bottom-right corner (black ellipse). 

Middle-right: CO(1-0) position- velocity diagram (PVD) taken over a one synthesized beam width (5 pixels) slice at the position angle 
indicated in the top-right corner inset and at the bottom of the panel (see £13.31 ). Velocities are heliocentric. The synthesized beam along the 
PV slice direction is indicated in the bottom-right corner. Contours are spaced by la starting at 3a (unless otherwise stated), where a is the 
rms noise in an individual channel multiplied by y/E, to account for the summation over the synthesized beam width. 

Bottom: CO(1-0) channel maps. The panels for velocity channels with flux are color-coded based on their velocity with respect to systemic 
(taken from Paper I). The mean velocity offset is listed in the top-left corner of each panel. The cross indicates a fixed reference position in 
all panels. Contours are spaced by la beginning at ±2cr (unless otherwise stated), where la is the rms noise in an individual channel. Grey 
contours are negative. The CARMA synthesized beam is shown in the lower-right corner of the first panel (black ellipse). 
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Figure Al. IC 676 is a field regular rotator (M^ = -22.27) with a bar and ring stellar morphology. It has a dust filament. The momentO peak is 17 Jy beam 
km s — 1 . The moment 1 contours are placed at 10km s _1 intervals and the PVD contours are placed at 1.5a intervals. 



20 K. Alatalo and the ATLAS 3D team 




Figure A2. IC 719 is a field 2a non-regular rotator (Mk = -22.70) with normal stellar morphology. It contains a dusty disc. The 2a peak is the signature of 
two counter-rotating stellar discs with opposite angular momenta. The molecular gas is co-rotating with a kinematically decoupled core, and counter-rotating 
with respect to the dominant stellar component. The momentO peak is 2.6 Jy beam -1 km s — 1 . The momentl contours are placed at 25km s — 1 intervals. 
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Figure A3. IC 1024 is a field regular rotator (Mk = -21.85) with stellar morphology that shows interaction features as well as a dust filament. The momentO 
peak is 14 Jy beam -1 km s -1 . The channel map and PVD contours are placed at 1.5cr intervals. 
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Figure A4. NGC 1222 is a field non-regular rotator {Mk — - 22.71) with stella r morphology that shows interaction features as well as dust filaments. It is a 
well-known interaction, and cross-identified as Markarian 603 (Beck et al. 2007). The momentO peak is 1 8 Jy beam 1 km s 1 . The PVD contours are placed 
at 1.5cr intervals. 
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Figure A5. NGC 1266 is a field regular rotator (Mk = -22.93) with normal stellar morphology. It is also observed to have a dust filament. Not only was it 
the brightest detection in the sample, but it had the largest linewidths observed. It was initially unresolved in CARMA D-array, thus was observe d in the much 
highe r resolution arrays. NGC 1266 hosts a massive molecular outflow that appears to be driven by an AGN. Further details can be found in lAlatalo et ail 
l201ll . The momentO peak is 100 Jy beam -1 . The momentl contours are placed at 10km s _1 intervals. The channel map contours are placed at 2a intervals 
and the PVD contours are placed at 3<r intervals. 
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Figure A6. NGC 2764 is a field regular rotator (M^ = -23.19) with stellar morphology that indicates interaction. It is also observed to contain a blue nucleus 
and dust filaments. The momentO peak is 18 Jy beam -1 km s _1 . The PVD contours are placed at 1.5a intervals. 
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Figure A8. NGC 3182 is a field regular rotator (M^ = -23.19) with normal stellar morphology and contains a dust bar and 
detections within the sample, thus it is likely much of the CO emission in this system is below the noise in the channel maps. 
beam -1 km s _1 . 



ring. It is one of the faintest 
The momentO peak is 3.3 Jy 
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Figure A9. NGC 3607 is a field regular rotator (Mk = -24.74) with normal stellar morphology. It contains a dust disc. The momentO peak is 17 Jy beam 1 
km s _1 . The moment 1 contours are placed at 40km s _1 intervals and the PVD contours are placed at 3<r intervals. 
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Figure A10. NGC 3619 is a field regular rotator (Mr- = -23.57) with a shell stellar morphology. It contains dust filaments, bars and rings. It is very likely that 
CARMA has resolved out what is possibly an extended gas disc. Observations that are sensitive to larger size scales is required to confirm this. The momentO 
peak is 6.9 Jy beam -1 km s -1 . The moment 1 contours are placed at 40km s _1 intervals. 
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Figure All. NGC 3626 is a group regular rotator (M^ = -23.30) that includes a double maximum velocity feature, with ring stellar morphology. It contains 
a dust disc. Unfortunately the central velocity of the observations were offset from Vsysi and due to its large linewidth, a small amount of the blueshifted 
spectrum was not within the 420 km s _1 CARMA window. The momentO peak is 4.3 Jy beam -1 km s -1 . The moment 1 contours are placed at 30km s -1 
intervals and PVD contours are placed at 1.5a intervals. 
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Figure A12. NGC 3665 is a group regular rotator (M^ = -24.92) with normal stellar morphology. It contains a dust disc. The momentO peak is 24 Jy beam 
km s — 1 . The moment 1 contours are placed at 50 km s — 1 intervals and the PVD contours are placed at 2a intervals. 
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Figure A13. NGC 4119 is a Virgo regular rotator (Mr- = -22.60) with normal stellar morpholog; 
km s — 1 . The PVD contours are placed at 1.5a intervals. 



;y. It contains a dust disc. The momentO peak is 8. 1 Jy beam 
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Figure A14. NGC 4150 is a Virgo regular rotator (M^ = -21.65) with normal stellar morphology. The momentO peak 
contours are placed at 1.5<r intervals. 



is 9.9 Jy beam -1 km s -1 . The PVD 
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Figure A15. NGC 4324 is a field regular rotator (M^ = -21.61) that includes two velocity maxima. It contains a ring stellar morphology, as well as a dust 
disc, bar, and ring. The momentO peak is 2.0 Jy beam -1 km s _1 . The PVD contours are placed at 2a intervals. 
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Figure A16. NGC 4429 

disc. The momentO peak 



is a Virgo regular rotator (Mk = -24.32) that includes two velocity maxima, with a bar and ring stellar morphology. It contains a dust 
is 14 Jy beam -1 km s _1 . The momentl contours are at 40 km s _1 intervals and the PVD contours are placed at 1.5a intervals. 
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Figure A17. NGC 4435 is a Virgo regular rotator (Mk = -23.82) that includes two velocity maxima, with a normal stellar morphology. It contains a dust disc. 
The momentO peak is 20 Jy beam -1 km s -1 . The moment 1 contours are at 30 km s _1 intervals. 
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Figure A18. NGC 4694 is a Virgo regular rotator (M^ = -22.15) with normal stellar morphology. It appears to be on its first approach into Virgo. It contains 
a dust filaments and bars. The momentO peak is 6.0 Jy beam -1 km s -1 . The moment 1 contours are placed at 10km s _1 intervals and the PVD contours are 
placed at 1.5a intervals. 
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Figure A19. NGC 4710 is a Virgo regular rotator (Mr- = -23.52) with normal stellar morphology. It appears to be on the outskirts of Virgo and is likely to be 
on its first approach. It contains a dust disc, and is one of the nearest to edge-on systems in the sample. The momentO peak is 62 Jy beam -1 km s -1 . Channel 
map and PVD contours are placed at 3a intervals. 
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Figure A20. NGC 4753 is a field regular rotator (Mk = -25.09) that includes two velocity maxima, with stellar morphology consistent with interaction. It 
contains a dust filament. It is also the most massive galaxy in the CARMA ATLAS 3D survey. The momentO peak is 12 Jy beam - 1 km s _ 1 . Moment 1 contours 
are placed at 40km s _1 intervals and PVD contours are placed at 1.5(7 intervals. 
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Figure A21. NGC 5173 is a group regular rotator (M^ = -22.88) with normal stellar morphology. It contains a dust bar. NGC 5173 was observed in CARMA 
D- and C- arrays, and the data were shared between this survey, and the work of Wei et al. 2010b, and thus uses the smaller pixel size of 0.4", to reflect the 
superior resolution of the observations. The momentO peak is 4.6 Jy beam -1 km s — 1 . PVD contours are placed at 1.5a intervals. 



40 K. Alatalo and the ATLAS 3D team 



45' 10' 



50' 



30' 



44' 10' 



59° 43 r 50' 



i 1 1 1 1 i 



NGC 5379 



13 h 55 m 37 s 








35 s 



33 s 



1500 1600 1700 1800 1900 2000 



0.3 
0.2 



3 
LL 



0.0 



■ 1 i 


■ 1 1 

fl WW 


\ 


K 


■ ■ i 

L 




7 1 


i 


up 

i i 


M 



■300 -200 -100 100 200 30C 
Velocity (km s" 1 ) 



59° 4410 




59° 4410 




13 h 55 m 35.5 s 33.5 s 



13 h 55 m 35.5 s 33.5 s 



10 20 30 40 50 
Offset (arcsec), p.a. = 66.° 



. -119 

■ ° © ■ 

+ 

' o ' 

: ■ 


-109 

o + 

o ^ 


-99 

c 

o + 

: ; © 


-89 

P o 

1 1 1 1 1 c> \ ' 


o -79 

' <3> ' 

: + ° <=; 

° 


-69 

■ O o - o- 

O i 


=> d : 59 ; 

+- 

° o ^ : 


' ° O 

+ 

; - o- 


g>-39 ■ 


■ <m -29°- 

■ ^ <m ■ 

o ° 


^ -19 

c? . 




1 


. - J| • : 

<r ; 


@>21 

o 

o 

G O ' 


^ ° - : 

■ D ■ 


■ o ® 

? - ; 

■ ^ " CD ■ 


: "°" o 51 : 

° ■ 


;<^ ft - 61 

■ ^ ' 

o 


71 

^ ° ^ o : 


:' ®- si 

;@ ® <*/ ; 

° n ° 

O °. 


91 

■» @ +° <a> - 

'■ (? : 

: ° s°: 


101 

o 

+ ° ■ 

o 




^ 111 

: ^ o : 

■ + o • < 

o 

: ;> - 8 : 


121 

+ 


^131 

' <^ ' « ° ' 

o 

o 

■ .0 . . _ 





40 r/ 
25 r/ 

59° 44' 10' 



13 h 55 m 35.5 S 33.5 S 



Figure A22. NGC 5379 is a group regular rotator (M K = -22.08) with ring stellar morphology. It contains a dust bar, ring, and filaments. The momentO peak 
is 5.7 Jy beam -1 km s -1 . PVD contours are placed at 1.5a intervals. 
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Figure A23. NGC 5866 is a field regular rotator (Mr- = -24.00) with normal stellar morphology. It contains a dust disc. The molecular gas also strongly 
suggests that this galaxy contains a bar. The momentO peak is 54 Jy beam -1 km s -1 . Moment 1 contours are placed at 40km s -1 intervals. Channel map 
contours are placed at 2<j intervals and PVD contours are placed at 2a intervals. 
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Figure A24. NGC 6014 is a field regular rotator (M^ = -22.99) with normal stellar morphology. It contains a dust bar, ring and filaments. The momentO peak 
is 14 Jy beam -1 km s -1 . PVD contours are placed at 1.5a intervals. 
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Figure A25. NGC 7465 is a field non-regular rotator (M^ = -22.82) with a kinematically-decoupled core (KDC) with stellar morphology indicative of 
interaction. It contains dust filaments. The momentO peak is 26 Jy beam -1 km s _1 . Channel map contours are placed at 2a intervals and PVD contours are 
placed at 4cr intervals. 
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Figure A26. PGC 029321 is a field regular rotator (Mk = -21.66) with normal stellar morphology. It contains dust filaments. The momentO peak is 8.5 Jy 
beam -1 km s — 1 . 
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Figure A28. UGC 05408 is a field regular rotator (M^ = -22.03) with a bar stellar morphology. It contains a dust bar and filaments. At a distance 45.8 Mpc, 
it is the most distant detection in the CARMA ATLAS 3D survey. The momentO peak is 3.5 Jy beam -1 km s _1 . 
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Figure A29. UGC 06176 is a field regular rotator (M^ = -22.66) with bar and ring stellar morphology. It contains a dust bar, ring, and filaments. The momentO 
peak is 14 Jy beam -1 km s -1 . 
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Figure A30. UGC 09519 is a field regular rotator (M^ = -21.98) with normal stellar morphology. It contains dust filaments. The momentO peak is 18 Jy 
beam -1 km s _1 . 
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Figur e Bl. Interfe rometric CO(l-O) data from ATLA S 3 D in the li terature. From top-t o-bottom, N GC 524 dCrocker et alJl201lh. NGC 2768 dCrocker et al 
2003) . NGC 3032 ([Young. Bureau & Cappellarill2008|). NGC 3489 JCrocker et all201lh . NGC 4459 jYoung. Bureau & Cappellarill2008l) . NGC 4476 iYoung 
2002I) . NGC 4477 dCrocker et alJl201lh . NGC 4526 JVoung. Bureau & Cappellarfcooi) and NGC 4550 JCrocker et alJl2009h ?For each galaxy, from left-to- 
right, the panels show the r-band image overlaid with the CO(1-0) integrated intensity (momentO) contours (cyan), the color-scale and contours of the CO 
integrated intensity (momentO) map overlaid with the IRAM 30m telescope beam (gray), the color-scale of the CO mean velocity (momentl) map overlaid 
with isovelocity contours and the IRAM 30m telescope beam (gray), and the g — r colour image overlaid with the CO(1-0) momentO contours (magenta). For 
comparisons to unsharp-masked HST data of these galaxies, see the papers cited above. 
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Figure Bl. Continued 
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Figure CI. NGC 2697 (Mk = -22.19) was removed from the sample as being mis-classified and containing spiral structure in the stellar isophotes. The 
momentO peak is 4. 1 Jy beam - 1 km s _ 1 . 
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Figure C2. NGC 4292 is a Virgo cluster galaxy (Mk = -23.2). It was removed from the ATLAS 3D sample based on its not being observed with SAURON, 
and thus lacks stellar kinematic data. The momentO peak is 4.3 Jy beam -1 km s _1 . 



